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Abstract 

In this study, internal energy (U), electric field (E) and particle 
number (N) which specify the system quantities i.e. thermodynamical 
quantities for the proteins. In the frame of thermodynamical formal- 
ism, the relation between the heat capacity at effective field Ce and 
the heat capacity at total dipole moment Cm and the relations for 
the increment of enthalpy AH, entropy AS and Gibbs energy AG 
which come out in the dissolving of the proteins in water have been 
obtained. By thinking about the present system being in the heat bath 
the canonical distribution function, by considering the system in heat 
and electric field bath the macro canonical distribution function and 
once more by taking the system in heat and particle number bath the 
macro canonical distribution function have been calculated. Partition 
functions have been related to the macro canonical quantities with the 
help of the free energy. Understanding the structure of the proteins 
have been endeavoured by fitting the theoretical calculations to the 
curves of experimental results. 
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1 Introduction 



Now a days, a progress in the investigation of the soft materials by the 
methods of statistical physics is recognized. It is not possible to think a life 
without the soft materials; all biological structures, molecules of the genetic 
code, proteins and membranes contain soft materials. The soft materials 
form the fundamental of life as well as being an important element of future 
technological developments. 

Statistical mechanical methods are the most practical methods in ob- 
taining the difference between the additional heat capacities which aries in 
the dissolution of proteins (one of the soft materials) in water and also for 
the calculation of thermodynamical functions. In favour of statistical ther- 
modynamics it is possible to explain the microscopic quantities by relating 
them to the macroscopic experimental results. 

While thermodynamical variables of fluids are represented by the phys- 
ical quantities U (energy), V (volume) and N (particle number), passing to 
the measurement results in the correspondence to T (temperature), P (pres- 
sure) and [i (chemical potential) consecutively In a similar manner, for 
the magnetic systems the thermodynamics variables are represented by the 
physical quantities U (energy), H (magnetic field) and N (particle number) 
passing to the measurement leads to a consecutive correspondence with T 
(temperature), M (magnetization) and n (chemical potential) 

For thermodynamical quantities of the proteins; internal energy (U), 
electric field (E) and particle number (N) could be taken as macroscopic 
variables. When the system tends to thermodynamics equilibrium, to relate 
with thermodynamics the system is determined by internal energy param- 
eter; temperature (T), electric field parameter; dipole moment (M) and 
particle number parameter; chemical potential (//). 

2 Relation between the heat capacities of proteins 

By taking electric field E and total dipole moment M as thermodynamics 
variables the first law of thermodynamics for the proteins could be written 
in the form: 



Here, the electric field is not an external field but is used to model the 
ice-like behavior of the water molecules around the non-polar surface. The 
electric field is a result of the effective behavior of the non-polar dissolvent 



dU = TdS - MdE. 
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that is applied to protein unfolding. By taking the proteins as a soft material, 
thermodynamical relation between the heat capacity C E at effective electric 
field and heat capacity Cm at constant total dipole moment has been put 
forward. Rearranging ^ one gets: 



dQ = dU + MdE. (2) 

To constitute such a relation if firstly the internal energy U is written as 
the total differentials of the variable T and E one could write 

U = U(T, E) dU(T, E) = (^)EdT + (^) T dE (3) 
and then substituting in Ogives 

dQ = (^)EdT + (^)rdE + MdE. (4) 

Rewriting 

dQ = (^) E dT+[M + (^) T ]dE. (5) 
When heat capacities are defined; since E = const, and dE = 0,0 becomes 

dQ = {^) E dT = C E dT. (6) 
Similarly when M = const, and dM = 00 then takes the form 

dQ = (^)EdT +[M + (^)r]dE = C M dT. (7) 

After substituting the definitions of heat capacities given by El and into El 
one obtains 

dU 

C M dT = C E dT +[M+ (g^^dE. (8) 

Dividing both sides of [HI by dT for M = const, and since = (§^)a/ one 
gets 

C M dT = C E dT + [M + Q) T }(^) M . (9) 

In order to express the term in the right hand side of this equation in terms 
of experimentally measurable quantities, [J could be written in the form: 
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1 M 

d S = -dU + T dE. (10) 

When the total differential of the internal energy U(T, E) is taken as given 
in El and substituted infTUl 

I dU dU M 

dS = -[(^) E dT + (^)TdE] + -dE (11) 

is obtained. Rearranging this equation 

dS = ^{ 3 ^) E dT + ^[(jj^rdE + f}dE. (12) 
Similarly the total differential of the entropy S(T, E) is written in the form 

S = S{T, E) =► dS(T, E) = ( W )EdT + (g^kdE. (13) 

Since corresponding terms in ^] and EH must be equal to each other: 

,&S\ l,dU. u dS, l„dU. 

<«f )E = T ( ffT^ and( -SE )T = T + M| ' (14) 

Writing down the property of equality of the second derivatives [A(^)e = 
^(|^)t] is substituted. ITU becomes: 

d A.dU. . d r l r ,dU. 

dE¥df )E] = dT { T [{ dE )T + M] - (15) 
Taking the derivatives inEl an d after simplification one obtains: 



T'dE^'df' E ~ ~ 7 T l{ -'dE' T ~ T 1 + TdT^dE' T + T^~dT^ E ' ^ ' 

Since this expression is a total differential it is independent of the order of 
the derivative. Rewriting 1161 

d dU 1 ( dU \ M d dU dM 

and then rearranging 

is obtained. Now substituting in |U] one finds: 
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C M = C E + T(^) E (^) M . (19) 

Using the change rule of the total differential (^^)e(^)m(§^)t = — 1 one 
could write (^)e = (§^)m(§^)t- Substituting this expression in H9l gives: 

C M = C E -T{ d -^) T [(^) M ]\ (20) 

By taking use of Maxwell relations written for soft materials (proteins) and 
with the purpose of expressing in terms of the measurable parameters from 
thermodynamics the following definition is given 



W T = { dT )M = aM (21) 



which transforms 1191 into: 



,dM ■ 



C M = C E -Tai I {—) T . (22) 



Similarly when the definition 



dE 7t 
is written down |2*H1 now takes the form 



C M = C E -T^. (24) 
7T 

3 Thermodynamics of enthalpy, entropy and gibbs 
energy increments 

The first law of thermodynamics for the soft materials has been written by 
n When S and M are taken as direct variables, change in enthalpy reads 

dH = TdS + EdM. (25) 

The enthalpy increment, after taking the total differentials of enthalpy H = 
H(S, M) and entropy S = (T, M), is written in the form: 

f) J-f f)J-f 
dH = { -dS hldS + { dM )sdM (26) 
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dS = Q) M dS + (—) s dM. (27) 
By substituting VM and [271 in 1231 

,dH. .dS. lrT1 r,dH, . dS . ,dH. , 

^ = < 55 )-^)^ T + t W T + W s]dM - (28) 

is obtained. The expressions T = (^)m, E = (§^)s and -jf- = (fjp)M are 
substituted in[2H]and than integration of the equation leads to the enthalpy 
increment: 



rT r M ( T) AC 

AH = / C M dT + / [£ + T{—) T ]dM. (29) 



lM(T t ) 

For the entropy increment; the total differential of entropy 5 = S(T, M) 
is taken and (^)m = ~{§§j)t and -jf- = (|^)m are substituted, which 
after integration gives: 

AS= [ T ^dT+ p T \§f)TdM. (30) 

•M T JM(T t ) OM 

By taking use of the first law of thermodynamics given in ^ which has 
been written for the soft materials and taking into consideration the param- 
eters U, E and N of the soft materials, change in the Gibbs energy could be 
obtained as 

dG = -SdT + EdM + fidN. (31) 

With a similar approach, the Gibbs energy increment is determined in the 
form 

rT fMrT) r N ( T ) 

AG = - SdT + / EdM + / fidN (32) 

JT t JM(T t ) JN(Tt) 

by taking the total differentials of Gibbs energy G = G(T, M) and total 
dipole moment M = M(T,E). These expressions together with E = (Jj|)t 

<dT> 



and -S = M are substituted inland then the equation is integrated. 



4 Establishment of the relation between thermo- 
dynamics and statistical mechanics 

Proteins is a common name of complex macromolecules which are formed by 
a collection of a great number of amino-acids and which play vital roles in 
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the life times of all of the living creatures. As a physical object, the protein 
molecule is a rather big polymer which is formed by thousands of atoms; in 
other words it is a soft matter such that from the physical point of view this 
is a macroscopic system. 

The three dimensional structure of the proteins depends on different en- 
vironmental factors inside the cell. Unfolding proteins return to their natural 
state when the environmental factors are removed, in manner of exhibiting 
the fact that the three dimensional structure of a protein solely depends on 
its information of the amino-acid regularity, the probable configurations of 
this three dimensional structure are in astronomically large order. In the 
thermodynamical limit (N — > oo, E —> oo and = const.) it is not pos- 
sible to make calculations by taking into consideration each particle in the 
macroscopic system formed from the proteins, even by using the computer 
existing on Earth with highest operation capacity. Thus combination of 
thermodynamics and statistical physics is a necessity. In the investigation 
of these type of system with the help of statistical physics, the existence 
of a great number of systems that are identical to the system under con- 
sideration is taken into account and using the statistical methods for these 
system, the most probable state in which the real system could enter is the 
task of determination. 

Statistical ensemble mean approach is one of the appealed methods in 
statistical mechanics. In order to constitute the statistical ensembles in 
the system that are interacting with its surroundings the system has to be 
taken into account together with its surroundings. Distribution functions are 
determined in the equilibrium state using the Lagrange multipliers method 
from statistical ensemble theory. 

4.1 A general approach to the calculation of distribution 
function, entropy and thermodynamical potential 

4.1.1 Distribution function and partition function 

For a general approach to the calculations, let the parameters U, E, N of 
the thermodynamical system under consideration, fluctuate together. In 
order words, the system could exchange energy, particles as well as electric 
field. Probability function P which gives the distribution of the system in to 
microstate is a function of n, E and N. Here n is a quantum number arising 
from the energy dependence. Number of particles N and electric field E also 
depend on energy. The entropy S of the system is given by: 
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S — — k P n ,E,N^ n Pn,E,N- 
n,E,N 



(33) 



The following equations for the total probability, the average energy U, 
average electric field E and average number of particles N of the system 
could respectively be written; 



p n,E,N = 1 



(34) 



n,E,N 



^ Pn,E,NUn t E,N = U 
n,E,N 



(35) 



^ Pn,E,NE n = E 



(36) 



n,E,N 



E Pn,E, N N n = N 



(37) 



Using the maximum entropy principle, the distribution function P u ,e,n 
of the system could be determined by Lagrange multipliers method. Thus 
the problem is; to make the entropy given bv 1331 maximum under the con- 
straints given by the 1341 1351 1361 and 1371 In the mathematical formalism of 
Lagrange multipliers methods, the expression; 

Q = ~k E Pn,E,NlnP n< E,N + Q!l(l - E Pn,E,N) + a 2 (U - ^ Pn,E,NU n ,E,N) + 
n,E,N n,E,N n,E,N 



a 3 {E- Pn,E, N E n ) + a 4 (N - £ P„ ijB ,iviV n ) (38) 




should be maximum, where a±, a 2 , 03, «4 are undetermined Lagrange mul- 
tipliers. Differentiating Q with respect to P n ,E,N and equating to zero one 
obtains: 



Pn,E,N = exp{-(a + (3Un : E,N + 0E n + jN n )} 



(39) 



where 




(40) 



Taking under consideration the condition in I34j IHHl leads to 
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expo = Z(P, 9, 7) (41) 
where Z(f3,0, r f) is the partition function of the system and is given by 

Z(/3,0 )7 ) = E exp-(PU n ,E,N + eE n + iN n ). (42) 

n,E,N 

Making use of \^ the distribution function now becomes 

^ exp-(J}U nEN + 6E n + 'yN n ) 

Pn,E,N= 2. , } (43) 

where (3, 8, 7 are functions of the parameters U, E and N which describe the 
system. In the rest of the manuscript the distribution function given bv 1431 
will be referred as macro canonical generalized distribution function. 

4.1.2 Entropy 

When the distribution function given by EH] is substituted in and taking 
into account that macro quantities U, E, N peak perfectly such that they 
could be taken instead of their average values, the entropy of the system 
becomes: 

S = kf3U + k9E + hr/N + klnZ(/3,e,-y) (44) 

4.1.3 Determination of the Lagrange multipliers 

With the purpose of obtaining the physical counterparts of the multipliers 
/?, 9 and 7, the definition of temperature T, total dipole moment M and 
chemical potential are made; 

If = (Qfj)E,N (45) 
M = T{—) U>N (46) 

" = - T{ 8N )U > E - (4?) 
Furthermore, after differentiation it is easy to show that; 



1 dZ 
Z^~dp 



^(—)e,-y = -U = -U (48) 
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Z v <9 7 )/3,7 ' ' 

In the above equations, instead of the average values U, E and N are taken. 
Substituting 1441 in 1451 and making use of EiSlone ends up with: 

»-w <51) 

Similarly substituting of liil in lifil and H71 respectively leads to: 

In the content of these calculations, distribution function 1431 takes form: 

exp-/3(U n E,N + ME n - /iN n ) , 
Pn ^ N = Z(T,M,ri (M) 

and for the partition function one obtains: 

Z n ,E,N = ^2 ex P~P( U n,E,N + ME n - fJ,N n ). (55) 

n,E,N 

On the other hand for the entropy given by 



„ U ME nN , , 

S=- + — -? T + klnZ (56) 



could be written. 



4.1.4 Thermodynamical potential 

Commencing with the entropy namely l5f)l and considering 

*(r,M,/x) = -kTlnZ (57) 
as thermodynamical potential one ends up with 

3>(T, M, n) = F + ME - fiN (58) 
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where F = U — TS is the Helmholzt free energy. could also be written 
in a different form: 

$(T,M,yu) = G-nN (59) 
where G = U — TS + ME is the Gibbs free energy. 

4.2 Obtaining the canonical distributions from the general- 
ized distribution 

The micro canonical, canonical and macro canonical distributions that are 
frequently met in the literature [3] can be obtained as the special cases of 
the generalized distribution, which is given bv 1431 

4.2.1 System in a heat bath: canonical distribution 

Thermodynamical systems which obey micro canonical distributions are 
systems isolated from energy, electric field and particle number exchanges. 
Therefore, energy has a unique value. 

In the thermodynamical systems governed by the canonical distribution 
there is energy exchange but the system is isolated from electric field and 
particle number exchanges. For this system a\ ^ 0, 02 7^ 0, 03 = 0, 
and «4 = 0. When the corresponding multipliers are substituted in 1431 the 
distribution function takes the form; 

Pn = ^exp(-pU n ) (60) 
and the partition function 1421 becomes 

Z = Y,exp-(3U n . (61) 

n 

One order hand the entropy expression given bv 1441 reduces to 

S = - + klnZ. (62) 

Meanwhile the thermodynamical potential of the system could be written 
as 

F = U -TS (63) 
where F = —kTlnZ is Helmholtz free energy. 
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Statistical mechanics plays the role of bridging the understanding of the 
physical quantities in terms of their microscopic counterparts. The energy 
of the proteins system is determined by: 

(64) 

The heat capacity at effective electric field has been calculated from the 
expressions of partition function and internal energy giving; 

8TI 

Ce = -WH-Qph- (65) 

Heat capacity at constant total dipole moment is obtained by making use 
of expression [I] 

C M = C E -T(^) T [(^) M ] 2 . (66) 

In a canonical ensemble, connection with thermodynamics is not estab- 
lished by entropy but with Helmholtz free energy. The following relations are 
written between Helmholtz free energy and the thermodynamical quantities; 

Entropy; 

dF 

S = -(-Qf)E (67) 

Total dipole moment; 

dF 

M = -(—) T . (68) 

4.2.2 System in a heat and electric field bath: macro canonical 
distribution 

Let us take into account the system where energy U and electric field E 
fluctuate but the number of particles is constant. For this case a,\ 7^ 0, 
«2 7^ 0, 03 ^ 0, and 04 = 0. In a similar manner to the above mentioned 
distributions, the distribution function 



Pn,E = lexp - (3(U n>N - ME n ) (69) 



the partition function 



E = J2^P-P(Un,N-ME n ) (70) 

n,E 



the entropy 
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S = ^ + ^ + klnZ (71) 
and thermodynamical potential 

$ = F + ME (72) 
are obtained, where <I> = —kTlnz*. 

4.2.3 System in a heat and particle bath: macro canonical dis- 
tribution 

In thermodynamical systems obeying the macro canonical distribution, the 
system is in a heat bath and material bath. For such a system oc\ ^ 0, 
«2 7^ 0, 03 = 0, and 04 7^ 0. When the corresponding multipliers are 
determined and substituting in 1431 distribution function takes the form: 

Pn,N = -Lexp - f3{U n)N - fJ,N n ). (73) 

For the systems obeying macro canonical distribution, 1421 which gives the 
partition function, reads 

H = J2 ex P ~ f3( u n,N - l*N n ) (74) 

n,N 

the entropy equation now becomes 

S=^-f + kl„Z (75) 
and thermodynamical potential of the system is 

$ = F - fiN. (76) 

In the macro canonical ensemble, in the first approximation, when the 
equation expressing the canonical ensemble is extended to macro canonical 
case without putting any restriction on the particle number the, following 
equation could be used jSJ; 

00 

*(T,i5, M ) = J2 W [exp(£-)Z(T,E,N = 1)] N . (77) 

As it is recognized, the partition function is summed over by weighing 
with Gibbs factor which has been used in the calculation of undistinguishable 
particles. When the sum in[J3is calculated; 
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$(T,E,z) =expzZ(T,E,N = 1) (78) 

is obtained. In this equation, /x being the chemical potential, the term 
z = exp{(3[i) corresponds to the pH of the system and it is a weighing 
factor. The energy of the system is given by 

U = -(^),, E (79) 

and it is determined by differentiation after substitution of the partition 
function. In terms of this grand potential, the expressions for the other 
thermodynamical quantities are given by; 
Entropy; 



Total dipole moment; 
Particle number; 



S = -(^),,E (80) 
M = -(§§ W (81) 

N = -(™)t,e. (82) 

Enthalpy, Gibbs energy and entropy increments and physical quantities 
e, b and fj, in H23 EH1 EZI an d EE1 that has been obtained in this semi phe- 
nomenological theory are determined by fitting experimental results given 
by 00. 



5 Conclusions 

In this study, for the soft matters in the frame of thermodynamics formalism, 
the relation between the heat capacities in 1241 AH enthalpy, AS entropy 
and Gibbs energy AG increments have been calculated and presented in 1291 
l3Ul and lS2*l respectivelv. Establishment of the relation between the thermody- 
namics and statistical mechanics ha been done with the help of the partition 
function and given by EH E21 E01 and EH1 Calculation of the partition func- 
tions for proteins (soft materials) with multi molecules exhibits difficulties. 
Therefore the partition function of a single molecules is obtained with the 
mean field approximation and then multiplied by N. The partition function 
is related to the parameters e, b and /i which determined the properties of 
the molecule. In El presented here; Cm is taken from experimental data, Ce 
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is theoretically calculated and the above mentioned parameters are obtained 
by fitting to the experimental results. 

On the other hand in the calculation of the increments of enthalpy AH, 
entropy AS and Gibbs energy AG, one proceeds in a similar manner. As an 
example, in the calculation of the entropy increments, as mentioned above 
Cm is taken from the experimental data and the term related with entropy 
is obtained from IH7I As a conclusion one could say that; the structure of 
the proteins (soft materials) might be understood by the thermodynamics, 
statistical mechanics and experimental studies triplet using this semi phe- 
nomenological approach. 

The significance of the electric field is to model the ice-like behavior ex- 
hibited by the water molecules around the non-polar surfaces which is a kind 
of stiffness jHj • Here this is not a real external electric field but the effective 
behavior of the non-polar dissolute enforced upon in the protein unfolding. 
It could also be possible to approach to this problem by considering the 
protein as a rigid body. However, in the proteins which are example of the 
soft matters, the benefit of determination of the expression which displays 
the thermodynamical relation between the heat capacity at effective elec- 
tric field Ce and the heat capacity at constant total dipole moment Cm is 
evident. 

As it could be seen from |^] which shows the thermodynamical relation 
between the heat capacity at effective electric field Ce and the heat capacity 
at constant total dipole moment Cm, the difference between the heat ca- 
pacities depends linearly on temperature. However results of experimental 
studies indicate a nonlinear dependence. This in turn implies the necessity 
of taking into account the internal structure of proteins. Ce given in 
could easily be calculated theoretically and if Cm expression is obtained 
experimentally then this enables the system to be investigated in terms of 
thermodynamical quantities. 
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